Topology and metastability in the lattice Skyrme model 
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Abstract 



We offer the Skyrme model on a lattice as an effective field theory — fully 

quantized — of baryon-meson interactions at temperatures below the chiral 

. phase transition. We define a local topological density that involves the vol- 

, umes of tetrahedra in the target space S"^ and we make use of Coxeter's 

' formula for the Schlafli function to implement it. This permits us to calcu- 

00 ' 

' late the mean-square radius of a skyrmion in the three-dimensional lattice 

' Skyrme model, which may be viewed as a Ginzburg-Landau effective theory 

for the full quantum theory at finite temperature. We find that, contrary to 
, expectations, the skyrmion shrinks as quantum and thermal fluctuations are 

enhanced. We ascribe this to a large number of metastable states that become 
Q ' accessible as the temperature is raised. 

> : 

. I. INTRODUCTION 

X 

5^ \ The Skyrme model is a popular model of the dynamics of pions and nucleons, incor- 

porating the former as its fundamental, pseudo-Goldstone field and the latter as topological 
solitons. The continuum theory has been widely studied via semiclassical techniques, giving 
a satisfactory phenomenology of low-momentum and low-temperature physics. In this paper 
we develop a lattice formulation of the model.Q 

Our motivation in employing a lattice cutoff is to overcome the limitations of the con- 
tinuum theory. In point of fact, the term "continuum theory" is misleading. The Skyrme 
Lagrangian is a non-linear sigma model with a four- derivative term, which makes it non- 
renormalizable in perturbation theory. This means that calculations of quantum effects must 
involve a short- distance cutoff. In working with chiral Lagrangians coupled to a local nucleon 



field, this is not a serious problem [B-|TO[ . One absorbs divergences into an ever-lengthening 



^For previous work see 
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list of counterterms; as long as external momenta are kept small, dimensional analysis limits 
the contributions of these counterterms to the results. When the nucleon is a soliton rather 
than a fundamental field, however, things are more difficult. The soliton's collective degrees 
of freedom are quantized separately, and the systematic development of higher-order quan- 



tum effects involves disentangling the pion field from these collective quanta [^. In any 
case, there is no way to limit the addition of higher-dimension terms to the Lagrangian as 
long as one must consider energy scales approaching the nucleon mass. 

We propose to turn these points to advantage, by considering the Skyrme model as an 
effective field theory. The cutoff, and the cutoff scheme (the lattice), are part and parcel of 
the specification of the theory. The free choice of terms in the Lagrangian is now a virtue, 
constrained only by phenomenological necessity. The procedure will be to pick a Lagrangian; 
to fix the lattice spacing; and finally to calculate any physical quantity of interest, going 
beyond perturbation theory and beyond semiclassical methods, since the lattice offers many 
more direct techniques. From a lattice point of view, the fixed cutoff is a virtue as well, 
since no continuum limit is necessary (nor is it possible). 

The most interesting feature of the Skyrme model is the stability conferred on the soliton 
by the topological conservation law and the four-derivative coupling. We wish in fact to 
isolate a single skyrmion in order to see how its properties change with temperature, density, 
etc. A recurring problem in lattice studies of solitons is the tendency of lattice dislocations 
to destroy topological stability or, in other words, to allow topology to slip through the 
lattice. We avoid this problem by choosing a lattice action that enforces continuity at short 
distances. In order to keep our numerical evolution within a single topological sector, we 
insist on a local updating scheme that should not nucleate smooth skyrmions. We find that 
these two ingredients suffice to stabilize the lattice skyrmion. 

The main technical development in this paper is the construction of an exact topological 
density that can be measured on a lattice field configuration. By "exact" we mean that the 
winding number, the sum of the density over the lattice, is always an integer. We cut the 
lattice into fundamental tetrahedra and map each tetrahedron into a curved tetrahedron 
in the target space. The winding number is then the sum of the signed volumes of the 
tetrahedra in S^. We calculate the volume of a spherical tetrahedron via a formula due to 



Coxeter |I2[, derived as a solution of differential equations first written down by Schlafii 

g. 



We summarize the continuum Skyrme model in Sec. II in order to establish notation. 
General considerations regarding continuity and topology on the lattice are presented in 
Sec. Ill, which concludes with formulae for the lattice action we employ. In Sec. IV we 
present our definition of the local topological density via Coxeter's formula for the volume 
of a quadrirect angular tetrahedron in S^. As a first application, we present in Sec. V the 
results of Monte Carlo simulation for the classical Skyrme model in three dimensions. Easier 
to simulate than the full four-dimensional path integral, this model may be regarded as a 
Ginzburg-Landau theory for the Skyrme model at finite temperature — a sort of effective 
theory for the effective theory. We use the topological density to calculate the mean-square 
radius of a skyrmion as a function of the couplings f3i and (32 of the lattice action. If we fix 
P2/P1, so that the form of the action is fixed, we find a multitude of metastable configurations 
of the skyrmion that are smaller than the ground state solution. As the "temperature" is 
raised by decreasing Pi and P2, these metastable states are made accessible to fiuctuations, 
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with the result that the skyrmion shrinks as it is "heated." Presumably the "temperature" 
of the 3d theory is an increasing function of the real temperature in the 4d theory that it 
approximates. Thus we reach the result that a quantized skyrmion shrinks as it is heated. 



II. CONTINUUM SKYRME MODEL 



The Euclidean action of the Skyrme model is 

f2 1 



s 



d X 



f 

16 1^1^ 32e2 



TT{[L„L,]f 



where the non-linear chiral field U is an SU{2) matrix, and we have defined 



= lU^d^U. 



{21) 



(2.2) 



We have omitted a mass term of the form Tr U and thus S is invariant under the SU{2) x 
SU (2) group of chiral rotations, 



U AUB\ with A,B e SU{2). 



(2.3) 



Classically, the symmetry is spontaneously broken. Finite-energy field configurations 
must tend to a constant at infinity, and one can use a symmetry rotation to make this 
constant the unit matrix, i.e.. 



f/(x) 



U as Ixl 



oo. 



(2.4) 



Using the Pauli matrices Tj, i = 1, 2, 3, we can write U in terms of new fields cr = (ctq, cTj) 
via 



t/ = (To + iaiTi. 



(2.5) 

U, we can identify cTj, 



where cr ■ cr = 1. By considering small fluctuations about U 
i = 1, 2, 3 with the Goldstone pion field. 

Equation ( |2.4] ) means that 3-dimensional space is compactified to the 3-sphere 5*^. Since 
U G SU{2) also takes values in S'^ [see Eq. ( p.5|) ], configurations t/(x) can be classified 
according to the homotopy group 7r3(S''^) = Z. Thus there is an integer winding number n 
that denotes how many times t/(x) covers the 3-sphere in field space as x is varied over its 
3-sphere. This winding number is topologically conserved, meaning that it cannot change 
under continuous deformation of the field t/(x). 

The n = sector includes the vacuum ^7(x) = U and perturbations around it. An 
example of an n = 1 configuration is the spherically symmetric skyrmion. 



U (x) = exp 



(2.6) 



with 
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/(O) = vr, 

/(oo) = 0. (2.7) 



/(r) should be determined so as to minimize the static energy 



E= 



P 1 

16 1*1^ 32e2 jj; 



(2.8) 



but a simple choice with the right topology is 



/(r) = TT (l -tanh— ) . (2.9) 
V ro/ 

Skyrme identified n with the baryon number of a field configuration, and the lowest soliton 
configuration with the nucleon. 

Given an arbitrary field configuration t/(x) satisfying Eq. ( p.4|) , its winding number may 
be calculated with the formula 

n = J dxeijkTr LiLjLk. (2.10) 

The geometric meaning of Eq. ( |2.10| ) will be apparent in its lattice counterpart below. 



III. LATTICE TOPOLOGY AND CONTINUITY 

In this section we define a topological density for a lattice field configuration that is 
unambiguous and that sums exactly to an integer. The further demand of conservation of 
the winding number will lead us to choosing a lattice action that constrains discontinuities 
in the field. 

A lattice configuration is specified by the field G SU (2) or, equivalently, by the 4- 
vector cTn G S^. In order to define the winding number, we begin ^ by cutting the cubic 
lattice into tetrahedra, five tetrahedra per cubic cell (see Fig. |l|). The four vertices n'^*\ 
i = 1-4 of each tetrahedron map to four unit 4- vectors cr(*\ which are vertices of a spherical 
tetrahedron in S^. This tetrahedron is defined via its vertices; its edges are arcs of great 
circles and its faces are spherical triangles drawn on great spheres. 

Since three (non-collinear) points in 5*^ determine a great sphere, it is clear that two 
adjacent tetrahedra in map into adjacent tetrahedra in S^, with the common face in 
mapping into a common face in 5*^. Thus the field configuration (Jn gives a triangulation 
of some volume in S^. We impose periodic boundary conditions on the lattice, and thus no 
tetrahedron possesses a face that is not shared by another tetrahedron (a free face). This 
implies that the complex of tetrahedra in 5*^ possesses no free faces either, and thus the 
volume covered by the complex must be an integer multiple of the volume^ of S^. If the 



^ The three-dimensional volume of the unit 3-sphere is 2tt'^. We will redefine this to be unity, and 
thus measure volumes in units of 2tt'^. 
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FIG. 1. A cube cut into 5 tetrahedra. The even vertices (circled) are connected to form a 
central tetrahedron and four others. In neighboring cubes, the odd vertices are to be connected 
instead, so that the diagonal edges match up. 

field configuration is smooth, the complex in will be composed of tetrahedra that are 
small compared to S^, and it will wrap around the sphere much as a smooth mapping of the 
3-torus must. In any case, the volume in S'^ of the five tetrahedra corresponding to a cube 
in gives a definition of the topological density pn contained in the cube. 

The topological density thus defined is not unambiguous, however. If a tetrahedron in 

is specified by its faces, then there are two volumes in that are bounded by these 
faces. One of the volumes includes the north pole (for example), and the other does not. If 
one of the volumes is measured to be a positive V (with V < 1), then the other volume will 
he V — 1. (Careful attention to the orientation will make the latter negative.) In order to 
assign a unique topological density to a lattice field configuration, we define a tetrahedron's 
volume to satisfy \V\ < |. 

We have solved the problem of uniqueness, but not that of conservation of the winding 
number. Consider a field configuration wherein one tetrahedron in has volume | — e, 
where e is small. Under a fluctuation of one of the vertices of the tetrahedron, its volume 
may shift to | + 5, while the neighboring tetrahedra change their volumes by — (e + 5) so 
that the winding number is unchanged. Unfortunately, our algorithm will now redefine the 
volume of the first tetrahedron to — | + 5, resulting in a loss of 1 in the winding number. 
This is often called "topology dropping through the lattice." 

The procedure to be followed at this point depends on the physics to be investigated. 
One might want to study, for example, the thermodynamics of the Skyrme model by fixing 
a chemical potential coupled to the winding number. The grand partition function would 
then be 

Z(/.)=^e'^"Z„, (3.1) 

n 

where Z„ is a sum over all field configurations with winding number n. A Monte Carlo 
simulation should then be allowed to wander freely among the sectors of different n, subject 
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to an acceptance/rejection test that enforces the relative probabihties e^". A local updating 
algorithm can nucleate "point skyrmions" as in the preceding paragraph, which change n by 
a unit and then spread out into smoother skyrmions; in addition, a non-local update could 
be permitted that creates smooth skyrmions directly. 

Our interest, however, is in the properties of a single skyrmion, which are accessible 
through a canonical ensemble at fixed n = 1. The simulation must be constrained so as 
not to change n. This requires that we prevent the nucleation of both point skyrmions and 
smooth skyrmions. The latter can be prevented by choosing a local updating algorithm. For 
the former, we choose a lattice action that excludes the possibility of tetrahedra in S'^ with 



1^1 ~ |. We adopt the kinetic term Wi 



Si = {a - log(crn ■ crn+^ - a) (3.2) 



in order to constrain cTn ■ crn+/i > a and thus to reject fluctuations that put large angles 
between neighboring field variables. With some Monte Carlo exploration, we can find a 
value of a that will keep the tetrahedral volumes far enough from |. 



For smooth configurations, we expand a-^ ■ o-n+ji around 1 and Eq. (|3.2|) becomes 

5*1 ~ ^(1 - (Tn ■ Crn+^) 

= E [l - ^Tr (u^Ui, 



= E (f/n+A - Un) (f/I+^ - Ul) , (3.3) 

nfj, ^ 

which approaches the kinetic term in Eq. ( |2.1| ) in the continuum limit. 

A technical point remains. We have discussed the ambiguity in fixing the volume of a 
tetrahedron if its faces are given, which led us to require \ V\ < |. A lattice field configuration 
CTn gives us only the vertices of each tetrahedron, however, not its faces. Given three vertices 
that determine a great sphere, the face that connects them can be chosen to be either of 
two triangles that together make up the sphere. A little thought shows that the difference 
between the volumes enclosed is |. This gives an ambiguity between a value V > and 
y — \ for the volume of the tetrahedron. 

This last ambiguity involves choosing which triangle constitutes the face of the tetra- 
hedron. But this face is shared between two adjacent tetrahedra, and the ambiguity will 
be immaterial if the same choice is made for both. This is hard to program, however, if 
each tetrahedron is to be handled on its own. We prefer to resolve the ambiguity for each 
tetrahedron separately, by requiring that \V\ < |. We impose this via a more restrictive 
choice of a in Eq. (|3.2|) . We find that setting a = 0.1 is adequate for the purpose. 

A local updating algorithm based on the action ( |3.2| ) will conserve winding number. Even 
so, a skyrmion will tend to shrink down to a point in accordance with the scaling argument 
that gives Derrick's Theorem. This is because the stabilization is being done at the scale of 
a single lattice spacing; only here will the scaling argument fail and the skyrmion run into a 
repulsive potential. In order to have a stable skyrmion of appreciable size, we cannot avoid 
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adding a Skyrme term as in Eq. ( |2.1| ). The most straightforward lattice transcription uses 
symmetric derivatives @], 

n ii>u 

(3.4) 

This couples each site to sites two links away, effectively via plaquette couplings across 
plaquettes of side \/2, and fails to couple the odd and even sublattices. We prefer to use a 
discretization that couples only across single plaquettes, of side 1: 



n fi>u 

(3.5) 

Our lattice action is a combination of the kinetic and Skyrme terms, 

S = piSi + P2S2. (3.6) 
In the naive continuum limit, we can compare with Eq. (|2.1| ) to deduce that 



IV. TOPOLOGICAL DENSITY ON THE LATTICE 

As shown in Fig. |l|, we cut each cube on the lattice into five tetrahedra. Each tetrahedron 
maps into a curved tetrahedron in 5*^, the volume of which gives the local topological density. 
We give here a practical, exact formula for this volume. 

In the 2-sphere 5"^, the area of a triangle is given simply by Girard's theorem, A = 
a + /3 + 7 — TT, where the right hand side is the sum of the triangle's angles minus vr, known 
as the angular excess. The easiest proof of the theorem is an argument based on overlapping 
lunes. This argument may be generalized [|T6| to a simplex in S*", but it only gives a formula 
for the volume when n is even; for odd n, one obtains instead a constraint on the angles of 
the simplex. 



L. Schlafli |13| attacked the general problem of the volume of convex polytopes in S'^' 



He was led to consider (for n = 3) the special case of the quadrirectangular tetrahedron, 



which we shall abbreviate as qrt. A qrt (see Fig. |^) is constructed as follows [0. Choose 
four points Pq, Pi, P2, P3 such that the line segments PqPi, P1P2, P2P3 are all mutually 
perpendicular. Then connect all four points together to form a tetrahedron. Every face of 
this tetrahedron is a right triangle. Three of the dihedral angles — the ones at P1P2, P1P3, 
and P0P2 — are right angles. The dihedral angles at P2P3, P0P3) and PqPi are not right 
angles, and they are denoted a, P, 7. 

This construction works equally well in Euclidean space and in S''^. A simple example of 
a spherical qrt is constructed by taking Pq, Pi, P2, P3 to be mutually perpendicular unit 
4- vectors. This is a qrt that covers 1/16 of 5*^. Its dihedral angles are all tt/2. 
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In Euclidean space, the dihedral angles of a qrt satisfy 

sin a sin 7 = cos/?. (4-1) 
In S^, they satisfy instead the inequality 

sin^ a sin^ 7 > cos^ /3. (4.2) 

The quantity 

D = \J sin^ a sin^ 7 — cos^ [3 (4.3) 

may thus be taken to be the generalization of the angular excess to this case. It vanishes as 
the qrt becomes small, which is the Euclidean limit. 

Let us denote the volume of a spherical qrt as V{a, /3, 7), normalized to 27r^ for the entire 
3-sphere. Schlafli derived formulas for the derivatives dV/da, dV/dp, dV/d'j (see also [|T7[] ). 
The integral of these formulas was found by Coxeter |]12[ in the form of a Fourier series valid 
in the restricted domain 

1 

< a< -TT 
- - 2 

< /5< TT 

< 7< -TT. (4.4) 
2 

Expressing the volume in terms of the complements of a and 7, 

y(a,/?,7) = ^5(|-«,/5,|-7), (4.5) 
his solution for the Schlafli function S is 
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, / D — sin X sin 2 \ COS 2mx — COS 27??,?/ + COS 2mz — 1 , , o , . 

S{x,y,z)=J2(-——. — 1 x^ + y^-z\ (A.Q) 

\D + sm X sm z / m"^ 

Armed with Eq. ( ^.61 ), we can calculate the volume of any tetrahedron in 5*^ by cutting 
it up into six qrt's (see Fig. |^). This is done by picking a vertex A of the tetrahedron and 



A 




C 

FIG. 3. A Euclidean tetrahedron cut into six qrt's. (Only two are shown.) 

dropping a perpendicular to the opposite face (the base) at E, then dropping perpendiculars 
from E to the edges of the base. Simple trigonometry gives the vertices of the qrt's in terms 
of those of the original tetrahedron; from the vertices of each qrt we calculate its dihedral 
angles a, /3, 7. Some technical problems, as well as further discussion of the Schlafii function, 
are relegated to the Appendix. 

The topological density pn in the lattice cube at n is the sum of the volumes of the 
tetrahedra in 5*^ that correspond to the Euclidean tetrahedra making up the cube. It goes 
without saying that summing pn over the lattice must give an integer (in units of 27r^) for 
any configuration, to high precision. 



V. CLASSICAL SKYRME MODEL IN 3 DIMENSIONS 

The minima of the action ( |3.6| ) in 3 dimensions are the static, classical solutions of the 
4-dimensional Skyrme model. When we include fluctuations about these minima, we can 
think of the 3d action that governs them as of an approximation to an effective action 
derived via dimensional reduction. Thus the parameters f3i, (32 are (unknown) functions 
of the parameters of the full 4d theory, namely, e, and the temperature T. Naturally, 
any reduction scheme that is precisely defined will give an effective action that is far more 
complex than ( |3.6| ), and further arguments will be necessary to justify its simplification. 

As an application of our algorithm for the topological density we calculate the mean- 
square radius of a single skyrmion in the three-dimensional theory, that is, of an equi- 
librium distribution of configurations with n = 1 at fixed f3i,f32- We fix initial conditions 
of the form ( |2.6| ) and (^.91), suitably discretized. Monte Carlo updates are done with the 
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usual local Metropolis algorithm, which preserves winding number as discussed above. The 
volume of the lattice is 16'^. 

The observable should be defined carefully. For each configuration used in the average 
one might calculate the barycenter of the topological density according to 



(5.1) 



and then the second moment via 



= ^min(n - Rc)Vn- 



(5.2) 



(The notation "min" means that one should take note of the periodic boundary conditions 
and always calculate the shortest distance between n and Rc.) Unfortunately, the boundary 
conditions make Rc calculated via Eq. (|5.1|) ill-defined. This may be illustrated by con- 
sidering a compact skyrmion with barycenter located on the face of the lattice at 2; = 0. 
Equation ( p.l|) fixes Rc in this case to be in the middle of the lattice, between the two 
half-skyrmions on opposite faces, and R^ will turn out to be on the order of the lattice size. 
The solution lies in regarding the lattice as a 3-torus on which the choice of origin (and of 
"faces") is entirely arbitrary. We evaluate Eqs. ( |5.1|) and (|5.2| ) for all choices of origin and 
take the minimal value of R^. The Rc that corresponds to this minimum is thus defined to 
be the location of the barycenter; this is as good a definition as any, considering the periodic 
boundary conditions. 

We show in Fig. the calculated skyrmion radius for f3i = 3 and various values of /32- The 
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FIG. 4. Mean-square radius of a single skyrmion, /3i = 3. The upper curve connects the radii 
of stable configurations; other points are radii of metastable skyrmion configurations. 



multiple values of R^ at a given value of /?2 refiect metastahility. These metastable skyrmion 
configurations are accessible from different initial states, specified by different values of tq 
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in Eq. ( |2.9| ). We have observed tunneling from metastable states to lower-energy states of 
larger radius that are apparently the global minima of the action. In all cases the lowest- 
energy states appear to be those of largest radius; these are shown connected by the upper 
curve in Fig. ^. The lower curve in the figure connects metastable states that were found to 
be mutually accessible by annealing. We have not explored the full set of metastable states 
for any given coupling; except for the annealing curve, the metastable states shown in Fig. ^ 
are effectively chosen at random by the initial conditions of the simulations. 

In Fig. ^ we present a more extensive set of results, this time showing variation with I3i 
as well as with 132- The curves connect data points for what we believe to be the equilibrium 
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FIG. 5. Mean-square radius of the topological density in the single-skyrmion sector 

skyrmion configuration at each coupling. Points not lying on curves are metastable config- 
urations of higher energy; again, we make no effort to show all such configurations at any 
coupling. We have not run across metastable configurations for /9i < 2. 

The abscissa in Fig. ^ is the ratio /32//3i that fixes the form of the action. The remain- 
ing overall coefficient j3i then acts as an inverse temperature that governs the fiuctuations. 
Regarding the 3d classical model as a Ginzburg-Landau theory for the fully quantized 4d 
theory, the lowering of (3i can be interpreted as enhancing the thermal and quantum fiuctu- 
ations in the latter. It is remarkable that heating the skyrmion by lowering I3i causes it to 
shrink. We ascribe this surprising property to the many metastable, excited states that are 
smaller than the ground-state configuration and that become accessible as the temperature 
is raised. The absence of metastability at the smallest values of (3i merely refiects the fact 
that the equilibrium skyrmion already averages over the various local minima. 

VI. DISCUSSION 

The qualitative conclusion of the preceding section must be taken with a large grain of 
salt. As mentioned above, we have no idea what the connection is between the couplings [3i 
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and of the 3d model on the one hand and the couphngs and temperature of the 4d theory 
on the other. It is possible that raising the physical temperature of the latter will lead to 
simultaneous changes in /3i and that cause the skyrmion to expand after all. In any 
case, the result is not quantitative and cannot be made quantitative until the dimensional 
reduction to the 3d theory is explored in detail. The metastable solutions, for that matter, 
might be artifacts of the specific lattice action we study that are absent in the true effective 
3d theory. 

Our study of the 3d model was motivated by its simplicity rather than by any funda- 
mental obstacle to studying the full 4d theory at finite temperature. Numerical work on the 
4d theory will merely require greater computer resources. In order to extract physics from 
lattice calculations, we will have to renormalize the theory by fixing the bare couplings (3i 
and (32, and the lattice spacing a, in terms of physical quantities such as the tttt scattering 
length and the skyrmion mass and radius. The temperature can be varied by changing 
the time extent Lf of the lattice in the time direction, or by varying the lattice spacing by 
changing (3i and (32 along lines of constant physics. The 4d quantum theory will offer as 
well the possibility of projecting the skyrmion to definite spin and isospin. Finally, it will be 
imperative to study the sensitivity of any physical quantities to the choice of lattice action. 
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We collect in this appendix some practical details concerning the use of the Schlafli 
function (|4.6| ) in calculating the topological density in lattice field configurations. 
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APPENDIX: 



1. Small qrt's 



Defining 



X = 



sin X sin z — D 



(Al) 



sin X sin z + D'' 



the formula for the Schlafli function is 



oo 



COS 2mx — cos 2my + cos 2mz — 1 



S{x,y,z)=Y.{-X) 



m 



711? 



2,2 2 

— a; + y — z . 



(A2) 



m=l 
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The sum converges absolutely when |X| < 1; since D > 0, this means sin x sin 2; > 0, which 
includes the entire domain (4^) considered by Coxeter. In the Euclidean limit, where the 
qrt is small, we have D ^ and hence X — 1. Then the sum in Eq. ( [A^ ) becomes the 
Fourier series for x"^ — y"^ + z^, and thus S" = in this limit as expected. 

The region of small D is where the series converges most slowly. We can avoid evaluating 
Eq. ( 1^) in this region as follows 0]. The vertices of a small tetrahedron in the 3-sphere 
are marked by 4-vectors cr^*) that are nearly equal. These 4-vectors are edges of a spindly 
4-pyramid with its apex at the center of the 4-ball bounded by the 3-sphere, the volume of 
which is 

n = l^e,.,.^H'a^MM''- (A3) 

The ratio of V4 to the total volume 7r^/2 of the 4-ball is approximately the ratio of the 
3-volume V of the tetrahedron to the total volume 27r^ of the 3-sphere, so 

V ^AVi. (A4) 

We calculate V4 for every tetrahedron obtained for each lattice cube (see Fig. 0). If it is 
small, we use Eq. ( |A4[ ) in lieu of cutting it into qrt's. We also calculate V4 for each qrt in 
turn, and if it small we use Eq. in lieu of evaluating the Schlafii function. Evaluation 
of V4 also gives an easy way to keep track of the sign of a qrt's volume, which is needed since 
S{x, y, z) is defined always to be positive. 



2. Large tetrahedra 

Another problem region for the Schlafii function arises when a ^ | or 7 — > |, giving 
X — s> — 1 and slow convergence. In the limit X = — 1, we have 

S'(x, z) = 7r(— X + y — z) = 7r(a; + /? + 7 — vr), (A5) 

a useful approximate formula in this regime. 

As we have noted, Coxeter's formula can only be used in the domain ([4.4|) . One can easily 
devise larger qrt's, however, and large qrt's arise regularly in calculating the topological 
density. We attack this problem by cutting such large qrt's into smaller tetrahedra. If either 
a or 7 is larger than ^, but not both, we determine the largest of 7 and draw a plane 
through its edge and the bisection point of the opposite edge (see Fig. |^). This creates 
two new, smaller tetrahedra. If both a and 7 are larger than |, we cut the tetrahedron 
by connecting all its edge centers; this gives four tetrahedra at the original vertices, plus a 
central octahedron. The octahedron can then be cut into four tetrahedra around any line 
connecting opposite vertices, of which we choose the shortest. 

Having cut the original qrt into two or eight pieces, we apply the algorithm (cut into 
qrt's and apply the Schlafii function) to each piece. Of course, any qrt that results might 
still be too large for Coxeter's formula. We then apply the cutting algorithm recursively. 
We have found that ten levels of recursion might be necessary for particularly rough lattice 
configurations. 
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Pi a P3 




FIG. 6. Cutting a qrt in the event that a is too large. 

The amount of recursion can be cut down dramatically by optimizing the operation of 
cutting the original tetrahedron into qrt's. There are four ways to choose the vertex A in 
Fig. 1^. For each choice of A, we determine the longest edge of the generated qrt's, measured 
by the largest pairwise angle among the 4-vectors cr*^*) of each qrt. This maximal edge is a 
measure of the size of the qrt's generated by the cutting process. We choose A in order to 
minimize the maximal edge. 
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